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PREFACE 

This  report  w.-is  prepared  by  Pieter  Hoekstra,  Geophysicist,  of  the  Physical 
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Engineering  Research  Branch,  Experimental  Engineering  Division,  and  Allan  J.  Delaney, 
Physical  Sciences  Technician,  of  the  Physical  Sciences  Branch,  Research  Division, 

U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory  (USA  CRREL). 

The  work  was  performed  under  DA  Project  4A162121A894,  Engineering  in  Cold 
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or  approval  of  the  use  of  such  commercial  products. 
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AIRBORNE  RESISTIVITY  MAPPING  OF  PERMAFROST 
NEAR  FAIRBANKS,  ALASKA 

by 

P.  Hoekstra,  P.V.  Sellmann  and  A.J.  Delaney 


INTRODUCTION 

Subsurface  investigations  for  geotechnical  projects  are  commonly  conducted  at  two  levels  of 
intensity: 

1.  Reconnaissance  investigations.  Information  is  often  obtained  from  existing  geological 
maps  and  airphotos,  which  provide  a general  impression  of  the  engineering  and  geological  aspects 
of  the  proposed  project, 

2.  Detailed  drilling  and  sampling  programs  at  selected  locations. 

In  the  mining  and  oil  industry,  it  is  customary  to  use  ground  or  airborne  geophysical  methods 
(seismic  and/or  electromagnetic)  after  reconnaissance  work  and  before  drilling.  Geophysical 
techniques  in  ,or  is  of  cost  and  in  terms  of  the  detail  of  information  they  provide  are  often  inter- 
mediate berwt  airphoto  mapping  and  sampling  by  drilling.  Ground  and  airborne  geophysical 
methtxls  attempt  to  derive  physical  properties  of  substrata  from  measurements  on  or  above  the  sur- 
face. Physical  properties  such  as  electrical  resistivity  and  seismic  velocity  often  correlate  with 
specific  soil  and  rock  types.  In  mineral  and  oil  exploration,  geophysical  techniques  are  often 
used  to  delineate  promising  targets  for  further  exploration  by  drilling  and  to  provide  a source  of  con- 
tinuous information  between  drill  holes.  Geophysical  techniques  now  only  find  limited  use  in  the 
gathering  of  subsurface  information  for  geotechnical  endeavors.  Two  reasons  for  this  are  apparent: 

1.  The  differences  between  the  physical  properties  such  as  electrical  resistivity  and  seismic 
velocity  in  near-surface  sediments  are  often  less  pronounced  than  the  differences  between  “cononnc 
mineral  deposits  and  surrounding  rocks. 

2.  A great  degree  of  sophistication  in  equipment  and  signal  processing  is  required  to  obtain  a 
continuous  profile  of  near-surface  physical  properties.  Therefore,  a contour  map  displaying  small 
nuances,  instead  of  identifying  anomalies,  is  desired. 

Electromagnetic  geophysical  techniques  are  the  subject  of  this  report.  The  development  of 
very  low  frequency  (VLF)  airborne  resistivity  mapping  techniques  by  the  U.S.  Geological  Survey 
(USGS)  (Frischknecht  1973).  Barringer  Research  Ltd.  (Barringer  and  McNeill  1969),  and  the 
Colorado  School  of  Mines  (Keller  et  al.  1970)  opened  the  possibility  of  obtaining  resistivity  data  at 
low  frequency  at  a relatively  modest  cost  and  in  inaccessible  areas.  The  latest  development  of 
multifrequency  resistivity  sensors  further  adds  to  the  potential  of  airborne  subsurface  mapping.  The 
North  American  Arctic  is  one  area  in  which  airborne  subsurface  mapping  can  fill  a large  need  for 
route  and  site  selection.  Therefore,  a study  was  conducted  with  the  Barringer  Research  Ltd. 
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multifrequency  airborne  E-phase  equipment  in  April  1973  in  the  vicinity  of  Fairbanks,  Alaska.  The 
objectives  of  this  study  were: 

1.  To  determine  the  ability  of  the  E-phase  system  to  delineate  various  surficial  geologic  units, 
in  particular  to  differentiate  between  perennially  frozen  and  thawed  ground. 

2.  To  investigate  in  detail  the  advantages  of  multifrequency  mapping,  which,  in  theory,  makes  it 
possible  to  examine  ground  conditions  at  various  depths. 

Several  study  sites  were  selected  in  the  vicinity  of  Fairbanks,  Alaska.  These  s4tes  were 
selected  for  several  reasons: 

1.  The  sites  are  in  the  discontinuous  permafrost  zone  and  consequently  have  a wide  range  of 
permafrost  conditions:  thawed  south-facing  hillsides  and  perennially  frozen  north-facing  hills, 
extensive  large  sections  of  thawed  ground  in  the  Tanana  Valley,  and  perennially  frozen  ground  in 
valleys  in  thf  ■'lands  near  Fairbanks.  The  permafrost  thickness  probably  varies  from  several  feet 
to  more  than  let. 

2.  The  geological  settings  and  associated  material  types  include'  a.  alluvial  silts,  sands 
and  gravels;  b.  colluvial  silts  and  organic  silts;  c.  eolian  silts;. and  d.  bedrock  that  ranges  from 
metamorphic  rock  to  intrusive  granitic  rock  types. 

3.  Substantial  amounts  of  background  information  on  the  local  geology  and  permafrost  distribu- 
tion are  available,  and  information  from  hundreds  of  logged  drill  holes  also  permits  local  and 
regional  correlations  to  be  made  (P^wd  1958,  Williams  et  al.  1 »59). 

The  location  and  the  layout  of  the  areas  surveyed  are  shown  in  Figure  1,  The  areas  surveyed 
with  the  E-phase  system  fall  into  two  categories:  l)  large  study  areas  that  cover  approximately 
20  square  miles  and  include  a range  of  surface  and  subsurface  conditions  commc  to  the  region, 
and  2)  small  sites  of  approximately  2 square  miles  chosen  because  of  the  av:  niity  of  detailed 

ground  control.  In  the  seven  areas  surveyed,  more  than  1000  line  miles  of  data  were  obtained  at 
two  or  three  different  frequencies.  In  this  report  the  results  of  four  areas  are  discussed:  Coldstream, 
site  2,  Chena  Hot  Springs  Road  and  Moose  Creek  Dam.  The  Rock  Quarry  Moose  Creek  Dam  is 
discussed  by  Sellmann  et  al.  (1974). 
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In  addition  to  the  airborne  survey,  ground  resistivity  measurements  were  made  at  selected 
sites,  and  laboratory  resistivity  measurements  were  made  on  selected  samples  to  help  translate  the 
resistivity  data  into  geological  formation. 

The  results  were  examined  on  a regional  scale  and  no  attempt  was  made  to  extract  quantitative 
information  on  local  variations,  such  as  permafrost  thickness,  and  gross  variations  in  ground  ice 
volume. 

Because  an  understanding  of  the  relation  between  resistivity  and  material  types  is  critical  to 
the  interpretation  of  resistivity  data,  resistivities  of  soils  and  their  variations  with  temperature 
and  other  factors  are  discussed  in  the  following  section.  In  the  section  Theory  and  method  (p.  7 ) 
the  pnticii  les  underlying  the  methods  are  discussed. 


RESISTIVITIES  OF  EARTH  MATERIALS 

Four  major  parameters  determine  the  resistivity  of  soil  strata'  (a)  soil  type,  (b)  watei  content, 
(c)  soil  temperature,  and  (d)  ice  content.  A large  number  of  less  important  factors.,  such  as 
pressure  and  the  free-salt  contents  of  the  pore  water,  also  influence  the  resistivity  of  soil  stiata. 
The  parameter  soil  type  is  independent  of  time  over  long  periods,  while  the  parameters  water  content 
and  temperature  are  subject  to  daily  and  seasonal  variations.  All  parameters  often  change  with 
depth  in  a profile  as  well  as  laterally. 


P (ohm-m) 


10'  I02  1O3  I04 


Figure  2.  The  resistivities  of  soil  types  in 
the  Unified  Soil  Classification  System  mea- 
sured in  exposures  in  the  valley  of  the  Tom- 
bigbee  River  in  Mississippi  (after  Hoekstra 
and  Delaney  1973). 


Dependence  of  resistivity  on  soil  type 

Figure  2 shows  a relation  between  resistivity  and 
soil  units  described  by  the  Unified  Soil  Classification 
System  (U.S.  Army  Engineer  Waterways  Experiment 
Station  1960).  The  measurements  weie  made  in 
exposmes  in  the  valley  of  the  Tombigbee  Rivei  in 
Mississippi.  The  clay  soils  (CH)  have  a lesistivity 
P less  than  50  ohm-m  and,  when  the  clay  content  of 
the  soil  decreases,  the  value  of  p uses  to  200  ohm-m 
in  typical  sandy  day  soils  (SC).  The  111-situ  amount 
of  water  in  the  exposures  probably  uinges  from  50 
to  100‘V,  saturation.  The  resistivity  of  giavels  also 
is  mainly  determined  by  the  amount  of  fines,  so  that 
clean  gravels  (GP)  show  resistivities  m excess  of 
1000  ohm-m,  while  gravels  containing  a few  percent 
oi  clay  (GC),  show  resistivities  as  low  as  200 
ohm-m.  The  relation  between  clay  content  and 
resistivity  is  influenced  by  the  dominant  type  of 
clay  mineral  in  the  soils,  and  the  actual  values 
within  a classification  change  from  one  area  to 
another. 


Relation  between  resistivity  and  water  content 

A typical  example  of  the  relation  between  water  content  and  resistivity  (the  inverse  of  conduc- 
tivity) is  shown  in  Figure  3 for  a o‘!t  soil;,  the  resistivity  p first  sharply  decreases  when  water  is 
added,  but  at  water  contents  in  excess  of  5%  the  rate  of  decrease  of  p with  water  content  levels  off. 
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Figure  3 . The  resistivity  of  a silt  soil  as  a function  of 
water  content  (g  H^O/g  soil). 

The  variation  in  p that  can  be  expected  in  situ,  because  of  naturally  occurring  changes  in 
water  content,  mainly  depends  on  the  groundwater  regime.  In  most  areas  of  Alaska,  the  water 
content  of  soil  probably  does  not  fall  below  the  water  content  at  field  capacity,  where  the  field 
capacity  is  defined  as  the  amount  of  water  held  in  soil  after  excess  water  has  drained  away  and 
the  rate  of  downward  movement  of  water  has  materially  decreased.  In  Fairbanks  silt,  resistivity 
changes  from  about  100  ohm-m  to  800  ohm-m  can  be  expected  because  of  changes  in  water  content. 

Dependence  of  resistivity  on  temperature 

The  relation  between  the  resistivity  of  saturated  soil  and  temperature  of  the  soil  is  given  in 
Figure  4 for  several  soil  types  and  one  rock  type  obtained  from  measurements  in  the  laboratory  by 
several  investigators  (e.g.  Ogilvy  1967,  Parkhomenko  1967).  At  temperatures  above  freezing, 
there  is  a small  increase  in  resistivity  with  decreasing  temperature,  and  it  is  only  below  freezing 
that  temperature  causes  substantial  changes  in  resistivity.  In  frozen  (round,  the  resistivity 
gradually  increases  when  the  temperature  is  lowered,  and  the  changes  ure  directly  related  to  the 
unfrozen  water  content  of  frozen  ground  (Hoekstra  and  McNeill  1973,  Anderson  and  Morgenstem 
1973).  The  general  trend  between  soil  type  and  resistivity  carries  over  in  the  frozen  state:  the 
clay-type  soils  have  the  lowest  resistivity,  followed  by  silt,  sand  and  gravel  soils.  The  resistivity 
of  frozen  soils  at  -5°C  is  about  a factor  of  5 higher  than  the  resistivity  of  the  same  soil  unfrozen; 
at  -1°C  no  more  than  a factor  of  1.5  separates  the  resistivity  of  frozen  and  thawed  clay  soils. 
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Figure  4,  The  resistivities  of  several  saturated  Figure  5,  Resistivities  of  saturated  organic  soils 

soil  types  and  one  rock  type  as  a function  of  of  Fairbanks,  Alaska,  as  a function  of  temperature, 

temperature. 


Soil  layers  rich  in  organic  content  are  common  in  the  first  5 m of  most  permafrost  ground,  and 
in  this  5 m of  surface  soil  the  annual  temperature  change  is  most  severe.  The  resistivity  on  some 
cores  from  surface  layers  of  permafrost  near  Fairbanks,  Alaska,  was  measured  as  a function  of 
temperature  in  the  laboratory.  The  relevant  data  on  these  saturated  soils  are  given  in  Table  1 and 
the  resistivity  as  a function  of  temperature  is  given  in  Figure  5.  Soils  high  in  organic  content  also 
often  have  high  water  content,  which  causes  a high  resistivity  when  soils  are  frozen.  Therefore, 
large  changes  in  resistivity  between  summer  and  winter  are  expected  for  the  organic  surface  soils 
in  permafrost  areas. 


Table  I.  Relevant  data  on  the  organic  samples  used  in 
measurements  of  data  shown  in  Figure  5. 


Sample 

Total 

density 

(g/cm*) 

Organic 

content 

(g/cm3, 

Water 
content 
( g/cm 3) 

Mineral 

content 

(g/cm3) 

Organic  silt  1 

1.02 

0.03 

0.82 

0.17 

2 

1.27 

0.04 

0.81 

0.42 

3 

1.37 

0.04 

0.91 

0.37 

Fairbanks  silt 

1.21 

<0.001 

0.61 

0.60 

Peat  1 

0.94 

0.12 

0.82 

0.003 

2 

0.94 

0.13 

0.95 

0.003 
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a.  Lenses  of  ice  up  to  7.5  cm  in  width. 


b.  Massive  structure  of  an  ice  wedge. 


Wi . 
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c.  Closely  spaced  small  lenses  of  ice  up  to 
2 mm  thick. 


d.  Pore  ice. 


Figure  6.  Different  ice  formations  in  the  USA  CRREL  tunnel  in  permafrost  in  Fox,  Alaska. 

Dependence  of  resistivity  on  ice  content 

Frozen  ground  in  permafrost  regions  can  have  ice  contents  far  in  excess  of  saturation,  and  r a 
lenses  and  massive  wedges  are  common  in  fine-grained  frozen  sediments.  The  USA  CRREL  tunnel 
in  Fox,  Alaska,  contains  exposures  of  silt  with  different  ice  types  and  contents  at  a temperature  of 
about  -3°C;  examples  of  ice  formations  in  the  USA  CRREL  tunnel  are  shown  in  Figure  6.  The 
relation  between  ice  content  and  resistivity  was  measured  by  a four-probe  method  with  the  probes 
inserted  in  the  tunnel  walls;  the  results  are  shown  in  Figure  7.  Frozen  silt  can  thus  have  a resis- 
tivity varying  from  about  1000  ohm-m  to  about  100,000  ohm-m  depending  on  the  ice  content. 


Resistivity  of  rocks 

The  resistivities  of  rocks  are  generally  more  difficult  to  classify  than  those  of  soils  because 
of  the  dependence  of  rock  resistivities  on  the  degree  of  fracturing  and  weathering,  grain  size  and 
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mineralogical  composition,  as  well  as  temperature  and 
water  content  of  the  rocks.  Some  in-situ  measurements  on 
rock  types  m the  Fairbanks  area  were  made  at  VLF  on  the 
ground.  The  Birch  Creek  schist  showed  resistivities 
ranging  from  100  to  2000  ohm-m,  and  the  local  intrusive 
rocks  displayed  resistivities  in  excess  of  2000  ohm-m 
(Sellmann  et  al.  1974).  Extensive  tables  on  rock  resistivi- 
ties, in  general,  are  given  by  Parkhomenko  (1967).  Heiland 
(1940)  and  Keller  et  al.  (1970). 


THEORY  AND  METHOD 

In  the  Barringer  E-phase  method,  the  electrical  prop- 
erties of  the  ground  are  derived  from  measurements  of  the 
groundwave  of  vertically  polarized  radiowaves.  This 
technique  has  been  known  for  a long  time  and  sevetal 
investigators  have  reported  on  its  use  (Eliassen  1956, 
King  1968,  Blomquist  1970,  Keller  et  al.  1970,  Hoekstra 
and  McNeill  1973). 

The  electromagnetic  field  vectors  of  the  groundwave 
at  grazing  incidence  are  shown  in  Figure  8.  There  is  at 
the  surface  a horizontal  electric  field  E^  in  the  direction 
of  propagation,  a horizontal  magnetic  field  Hy  perpendicu- 
lar to  the  direction  of  propagation,  and  a vertical  electric 


Figure  7.  The  relation  between  ice 
content  and  Resistivity  measured  in- 
situ  in  the  walls  of  the  USA  CRREL 
tunnel  in  Fox,  Alaska. 

ground  surface. 

agates  vertically  into  the  ground  with  field  vectors  Ex  and  Hy.  Beneath  the  surface  of  a homoge 
neous  ground  the  horizontal  electric  field  components  decay  with  increasing  depth  z as 


field  Ez.  BothE 


x and  Hy  are  continuous  across  the 
Thus,  at  gt  .zing  incidence  a wave  prop- 


Ex(z)  - E*(0)  e'V2 


(1) 


where  y is  the  propagation  constant  in  the  earth  given  by  (e.g.  Kraichman  1970): 


where  / - (~l)l/2 

p0  = magnetic  permeability  of  free  space  - An  * 10"  7 henry/m 
(0  - permittivity  of  free  space  = 1/36"  x 109  farad/m 
oj  - angular  frequency  (rad/s) 
p = ground  resistivity  (ohm-m) 
tr  = ground  relative  dielectric  constant. 

The  terms  j (fi0<o/p)  and  (<re0/i0<u2)  in  eq  2 are  the  conduction  currents  and  displacement  currents, 
respectively,  and  in  many  situations  encountered  in  the  ground  and  at  frequencies  below  10®  Hz 
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R«mot* 

Tron*mitt*f 


Figure  8.  Electromagnetic  Held  components  of  a vertically  polarized  radio  sur- 
face wave. 


Figure  9.  Skin  depth  of  electromagnetic  plane  wave?  as  a function  of 
frequency  and  ground  resistivity. 


P 


» i VoUo*'2 1 • 


(3) 


The  skin  depth  of  the  wave  propagating  into  the  ground  is  defined  as  the  depth  in  which  the  electro- 
magnetic field  decays  to  e-1  of  its  value  at  the  surface.  The  skin  depth  5 is  given  by 

S = Re(yfU/jt_Y'/\  (4) 

\m0W 


where  Re  (y)  indicates  the  real  part  of  the  propagation  constant. 
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Th1  skin  depth  is  a significant  parameter  in  that  it  indicates  the  effective  depth  of  electromagnetic 
wave  penetration  into  the  earth,  and  thus  the  effective  depth  to  which  resistivity  is  sensed,  since 
large  changes  in  resistivity  occurring  at  depths  greater  than  1 skin  depth  do  not  greatly  affect 
the  field  vectors  at  the  earth's  suiface.  Figure  9,  which  shows  the  behavior  of  the  skin  depth  as 
a function  of  frequency  and  ground  resistivity,  illustrates  that  the  effective  depth  of  investigation 
can  be  altered  by  changing  the  frequency.  Conversely,  a survey  may  be  made  with  several  fre- 
quencies, simultaneously  sensing  the  resistivity  to  selected  depths. 

Local  changes  in  the  electrical  properties  of  the  ground  cause  perturbations  in  amplitude  and 
phase  of  the  horizontal  electric  field.  Measurements  of  these  perturbations  at  or  near  the  surface 
are  the  basis  of  the  method.  Perturbations  in  Ez  and  Hy  occur  in  the  vicinity  of  fairly  abrupt 
lateral  changes  in  the  conductivity  of  the  earth,  but  these  fields  over  a limited  area  provide  no  in- 
formation about  the  properties  of  the  earth  when  gradual  variations  take  place.  However,  the  hori- 
zontal electric  field  is  dependent  on  the  electrical  properties  of  a horizontally  stratified  ground  and 
is  useful  in  mapping  gradual  variations  in  resistivity  and  thickness  of  flat-lying  strata  (Frischknecht 
1973). 

Diurnal  and  temporal  fluctuations  in  field  strength  occur  in  accordance  with  events  in  the 
ionosphere.  To  eliminate  the  effects  of  temporal  variations  in  field  strength,  the  horizontal  electric 
field  Ex  is  normalized  by  the  orthogonal  component  of  the  magnetic  field  Hy  to  obtain  the  surface 
impedance  Zs 


or  by  the  vertical  electric  field  £z  to  obtain  the  wave  tilt  W 


(5) 


(6) 


This  normalization  by  Hy  or  Ez  works  for  resistivity  mapping  because  Ez  and  Hy  are  virtually 
unaffected  by  gradual  resistivity  changes.  Both  Z8  and  W are  characteristic  of  the  refracted  wave 
penetrating  into  the  earth  at  the  location  of  the  receiver  and  are  independent  of  the  resistivity 
variations  found  in  the  path  traveled  by  the  wave  from  transmitter  to  receiver. 

The  theory  which  relates  the  measurement  of  the  surface  impedance  and  the  wave  tilt  to  the 
resistivity  of  homogeneous  or  layered  ground  has  been  well  developed  (e.g.  Wait  1962.  Kraichman 
1970).  For  a homogeneous  half  space  Zs  is  given  hy 


where  q - (1/p  + ;Vr<0<u)]'\  the  characteristic  impedance  of  ground 

y0  - propagation  constant  of  free  space  (p0<0)1* 
p - resistivity  of  ground  (cbm-m) 

9 = angle  of  incidence. 

For  all  values  of  frequency  and  dielectric  constants  likely  to  occur 
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Y0 


so  that  eq  7.  at  grazing  incidence,  reduces  to 


li- 


ft) 


The  relation  between  the  wave  tilt  W and  the 
surface  impedance  is  given  by 


V 


V*  y 
E/Hy 


l 

% W/ 


(9) 


where  the  characteristic  impedance  of  free 
space,  has  been  substituted  for  the  ratio  e /H 

T J 

In  many  situations,  conductive  currents 
dominate  displacement  currents,  so  that  eq  9 
can  be  approximated  by 


-13' 


UpHqoj)  ' (1+  /) 


Pff 


(10) 


Figure  10.  Magnitude  of  wave  tilt  as  a function 
of  frequency  for  resistivity  variations  over  typi- 
cal homogeneous  earth.  The  wave  tilt  has,  thus,  equal  in-phase  and 

quadrature-phase  components  for  homogeneous 
ground  if  displacement  currents  can  be  neglected.  This  also  means  that  the  phase  of  Ex  leads  that 
of  Ez  by  45°.  Figure  io  gives  the  wave  tilt  in  percentage  as  a function  of  frequency  for  resistivity 
variations  encountered  on  earth.  The  wave  tilt  is  in  most  instances  less  than  io%.  so  that  Ez  is 
often  at  least  an  order  of  magnitude  larger  than  F.  . 


In  the  E-phase  equipment,  Ex  and  Ez  are  measured  by  two  perpendicular  short  dipole  antennas. 
In  an  aircraft,  roll  of  the  plane  frequently  causes  leakage  of  £z  into  the  horizontal  antenna,  but  all 
of  the  leaked  signal  is  in  phase.  Therefore,  by  measuring  the  quadrature-phase  component  of  W and 
basing  the  computation  of  resistivity  only  on  that  component,  errors  caused  by  leakage  are  avoided. 
Hence,  resistivity  p can  be  computed  from: 


P 


<0w 


(ID 


where  V d indicates  the  quadrature-phase  component  of  i*'. 

Many  instances  are  encountered  m which  the  in-phase  and  quadrature-phase  components  of  the 
wave  tilt  are  unequal.  In  these  instances  an  apparent  resistivity  pa  is  defined  by  eq  u.  It  is 
called  an  apparent  resistivity  because  when  the  in-phase  and  quadrature-phase  components  of  the 
wave  tilt  differ  pa  is  unequal  to  the  true  resistivity  of  the  ground.  Often  the  true  resistivity  cun 
be  arrived  at  by  computation. 

In  permafrost  situations,  high  resistivity  values  can  be  encountered,  and  neglecting  displace- 
ment currents  can  result  in  errors  in  the  computation  of  pa;  this  error  is,  moreover,  amplified  by 
using  the  quadrature-phase  component  only.  Figure  li  shows  the  ratio  of  the  true  value  of  resistiv- 
ity  prea,  and  pa  as  a function  of  frequency  for  several  values  of  preai  for  homogeneous  ground;  a 
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Figure  It.  The  ratio  of  apparent  resistivity,  computed  from  eg  u 
fur  homogeneous  ground,  and  true  resistivity. 


dielectric  constant  of  20  was  assumed  in  all  situations.  For  homogeneous  ground  with  a resisti- 
vity of  1000  chm-m,  there  is  a close  correspondence  between  pa  and  preal  at  VLF.  but  at  broadcast 
band  (BCB)  frequency  fja  can  exceed  preal  by  as  much  as  a factor  1.3.  At  higher  resistivities,  the 
errors  become  progressively  worse,  and  at  a resistivity  of  105  ohm-m,  the  E-phase  yields  at  200  kHz 
a resistivity  of  about  101  ohm-m.  Fortunately,  uniform  ground  with  resistivity  values  in  excess  of 
IQ4  ohm-m  is  not  common,  and  often  there  is  little  practical  use  for  differentiating  between  resistiv- 
ities in  excess  of  5000  ohm-m. 

When  the  earth  is  not  homogeneous,  but  layered,  and  consists  of  several  layers  of  different 
resistivities  and  thicknesses,  the  surface  impedance  can  be  written  as  the  product  of  the  charac- 
teristic impedance  of  the  surface  layei  1 /t  and  a correction  factor  Q accounting  for  layering  (e.g. 
Kraichman  1970). 


Z vrQ. 

For  example,  for  a ground  that  is  stratified  in  three  layers,  Q takes  the  form' 

q _ )',0»  y.tanh^d, 

Y2  ~ y j^tanh  y,dj 


(12) 


where 


Q = 


y2  - y3  tanh  V2da 

y3  + ya  tanh  y2d2 


(13) 


and  y j , y2,  and  y3  are  the  propagation  constants  of  the  tliree  layers;  d j and  d are  the  thicknesses 
of  media  1 and  2;  and  the  thickness  of  medium  3 is  considered  of  infinite  depth. 
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Frequency  (Hz) 


a.  p1/p2  = 0.I. 


Frequency  (Hz) 

5.  Pj/pg  = 10.0. 

Figure  12.  The  computed  in-phase  / and  quadrature- 
phase  Q components  of  the  wave  tilt  and  its  phase  <t> 
as  a function  of  frequency  for  the  layered  ground  condi- 
tions shown. 

The  wave  tilt  over  layered  ground  has  unequal  in-phase  and  quadrature-phase  components;  in 
Figure  12  values  of  the  in-phase  / and  quadrature-phase  Q components  of  W and  its  phase  4>  are 
given  for  typical  two-layer  situations.  In  one  case,  the  ratio  of  the  resistivity  in  media  1 and  2, 
Pj/p2,  is  0.1,  typical  of  sites  where  soil  overlies  a granitic  bedrock,  or  where  silt  of  clay  over- 
lies  a gravel  deposit;  in  the  other  case,  pj/p2  is  10,  which  is,  for  example,  found  in  permafrost 
situations. 

The  phase  of  W is  quite  diagnostic;  when  px>  pz,  the  phase  is  greater  than  45°;  and  as  a 
corollary,  when  p}  < p2,  <f>  is  less  than  45°.  Unfortunately,  the  airborne  E-phase  system  does  not 
measure  the  phase  of  w. 
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Figure  13,  The  ratio  ot  computed  apparent  resistivity 
pa  and  the  resistivity  ot  the  surface  layer  pl  as  a func- 
tion of  thickness  ot  first  layer  in  terms  of  skin  depth. 


Figure  14.  Coverage  for  well  known  VLF  stations.  Circles  enclose 
approximate  areas  where  field  strengths  exceed  100  pV/m. 


Call 

Location 

Frequency  IkHti 

NAA 

Cutler,  Maine 

17.8 

NLK 

Jim  Creek,  Washington 

18.8 

NSS 

Annapolis,  Maryland 

81.4 

NBA 

Balboa,  Canal  Zone,  Panama 

84.0 

QBR 

Rugby,  England 

16.0 

UMS 

Corky,  USSR 

17.15 

NWC 

Northwest  Cape,  Australia 

83.3  and  15.5 

NPM 

Lualualei,  Hawaii 

83.4 

\ The  apparent  resistivity  of  a two-layer  ground  measured  with  the  system  depends  on  the  fre- 

quency; in  particular,  it  depends  n the  ratio  of  the  thickness  of  the  first  layer  and  the  skin  depth 
1 of  the  radiation  in  the  first  layer.  The  apparent  resistivities  of  two-layer  situations  are  shown  as 

a function  of  the  ratio  d xl  8 in  Figure  13.  In  Figure  13  the  apparent  tesistivities  are  computed  in 
two  ways;  p^^j,  based  on  the  quadrature  phase  component  of  W;  and  pabs,  based  on  1^1  • The 
value  of  pt  is  generally  bounded  by  the  resistivities  of  the  two  layers,  although  pflUtd  falls  at  cer- 
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tain  frequencies  below  the  value  of  Pl  for  pl  /pg  = 0.1.  The  trends  in  the  values  of  p uad  and,pat)8 
are  similar:  PAbs»p,uad  for  p j/pg  = 0.1  and  the  opposite  occurs  for  Pj/p2  = 10.  QUa 

Since  both  Pqy^jj  and  pubs  are  apparent  resistivities,  and  as  long  as  it  is  understood  that  an 
interpretation  requires  fitting  a solution  to  the  experimental  data,  the  behavior  of  the  curves  in 
Figure  13  shows  that  pflUtd  is  a satisfactory  parameter  with  which  to  characterize  the  resistivity 
of  a layered  medium.  A multifrequency  survey  yields  2 or  3 experimental  points  and  a satisfactory 
interpretation  of  the  layering  in  the  gjound  can  often  be  made.  If,  moreover,  through  other  inde- 
pendent measurements  the  values  of  pl  and  pg  are  known,  a realistic  interpretation  of  the  ground 
profile  can  often  be  made. 

The  existence  of  a strong  theoretical  base  for  the  E-phase  method  is  also  used  in  preliminary 
work  to  decide  on  the  merits  of  an  airborne  survey.  A general  idea  ot  the  range  of  resistivities  and 
geologic  conditions  in  an  area  can  be  obtained  on  the  ground  in  a matter  of  a few  days,  and  using 
these. resistivities  in  computer  modeling,  an  evaluation  can  be  made  of  the  ability  of  the  E-phase 
method  to  meet  certain  geotechnical  objectives.  A typical  example  of  such  an  evaluation  is  the 
work  in  the  area  of  the  Tennessee-Tombigbee  River  (Hoekstra  and  Delaney  1973). 

Surface  waves  of  the  type  described  above  are  generated  by  vertical  dipole  antennas  situated 
on  the  earth’s  surface.  Powerful  VLF  transmitters  situated  around  the  world  provide  sufficient 
field  strength  for  wave  tilt  measurements  in  the  range  of  15-25  kHz  over  a large  part  of  the  land 
area  of  the  world.  Figure  14  shows  coverage  for  well  known  VLF  stations.  The  circles  around 
each  station  indicate  the  areas  in  which  there  is  sufficient  field  strength  for  measurements.  In 
Alaska  the  station  NLK  operating  at  18.6  kHz.  and  situated  in  Jim  Creek,  Washington,  was  used. 
Military  transmitters  are  often  found  within  the  frequency  band  50  to  100  kHz,  and  navigational 
transmitters  are  often  found  within  the  frequency  band  200  to  400  kHz.  Thus,  in  most  instances 
surveys  can  be  made  at  three  frequencies  ideally  separated  for  effective  resistivity  sounding. 
Moreover,  in  the  band  200  to  400  kHz,  inexpensive  local  transmitters  can  be  erected  near  the  sur- 
vey site,  if  necessary. 


THE  E-PHASE  SYSTEM 

A horizontal  and  a vertical  antenna  are  installed  in  a cone  extending  from  the  nose  of  a short 
takeoff  and  landing  (STOL)  aircraft  (Britton  Norman  Islander)  as  shown  in  Figure  15.  The  E-phase 


Figure  15.  STOL  aircraft  (Britton  Norman  Islander)  with  E-phase  antennas  in  the  nose  cone  and  a 

magnetometer  stinger  on  the  tail. 
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instrument  is  capable  at  measuring  the  vertical  and  horizontal  components  of  the  electric  field  in 
four  frequency  ranges.  Of  the  four  ranges,  only  three  are  used  simultaneously  on  any  one  survey. 

The  instrument  can  be  continuously  tuned  within  the  frequency  bands,  and  phase-lock  loop  cir- 
cuitry is  utilized  to  keep  the  receiver  on  the  selected  frequency.  The  frequency  bands  are  as  follows: 

VLF  - range:  15*25  kHz 

SLF  * range:  70-140  kHz 
LF-  range:  200-400  kHz 
BCB*  range:  550*1100  kHz 

After  amplification,  the  signals  are  recorded  on  an  analogue  chart  and  on  magnetic  tape.  The  fol- 
lowing components  can  be  measured  at  each  frequency: 

Ez  (/)  - In-phase  component  at  the  vertical  electric  field. 

Ex  (/)  = In-phase  component  of  the  horizontal  electric  field. 

Es  (Q)  = Quadrature-phase  component  of  the  horizontal  electric  field. 


Calibratioa 


Calibration  of  the  EX(Q)  and  Ex(l)  channel  is  accomplished  by  tilting  the  horizontal  antenna. 

The  EX(Q)  channel  is  switched  to  read  Ex(l)  and  the  antenna 
is  tilted  20°.  The  signal  picked  up  in  the  tilted  antenna  is  Ez 
sin  20°,  and  the  percentage  of  wave  tilt  Ik  it  represents  is:. 


Olfaction  of  Flight 


E.-VLF 


Animator 


w 


E.  sin  20° 


x 100  = 34%. 


E,<QI-VL* 


E,(l)-tF(C00r»0) 

El(l)*BCB{Fin«) 


- -—Manual  Fiducial* 


E.-LF 

E.-BCB 


The  Ex(l)  and  EX(Q)  channel  RF  gains  are  set  during  calibra- 
tion to  give  a reasonable  deflection  on  the  chart  recorder  for 
the  20°  rotation  of  the  antenna,  and  changes  in  output  gain 
settings  during  flight  are  recorded  on  the  analogue  charts. 

The  calibration  of  the  signal  strength  of  the  horizontal  antenna 
in  terms  of  wave  tilt  is  dependent  on  the  signal  strength  of 
Ez(/).  If  Ez(l)  changes  during  the  course  of  the  flight,  then 
a correction  must  be  made  accordingly.  The  output  gain  at 
the  Ex(l)  channel  during  calibration  is  directly  proportional 


o' 


y'»  Agtomotie  Camara  Fiducial* 


Figure  16.  Staple  record  of  the 
analogue  chan  used  in  flight  re- 
corder. The  multiplexing  of  the 
Et  and  Ex  channels  is  shown  in 
coarse  and  fine  traces. 


Analogue  recorder 

The  signals  are  recorded  on  a six-channel  recorder.  A 
built-in  multiplexer  has  the  capability  in  this  system  of  multi- 
plexing all  6 of  the  chart  recorder  channels,  thus  doubling  its 
recording  capability.  This  multiplexer  has  been  set  to  switch 
at  approximately  twice  per  second  with  a mark  space  ratio  of 
approximately  5 to  1.  This  gives  two  very  distinct  traces  and 
removes  a great  deal  of  confusion  (traces  are  defined  as  fine 
or  coarse).  A sample  record  of  the  analogue  chart  (Fig.  16) 
indicates  the  designation  of  each  trace. 

Magnetic  recorder 

The  data  generated  during  a survey  flight  are  also  recorded 
on  a 24-channel  magnetic  tape  recorder  manufactured  by 
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Table  II.  Chauel  designations  of  Magnetic  tape  recorder. 


Ttpe 

recorder 

channel  Parameters 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

19 

13 

14 
16 
16 

17 

18 

19 

20 
21 
22 

23 

24 


Real-time  clock  - TH-UH-TM* 


Real-time  clock  • UM-TS-US* 


Bz(l)  - VLF 
Bx(l)  - VLr 
BX(,Q)  • VLF 
Hy(l)  • VLF  - 
Hy(l ) - VLF  - 
Hy{Q)  - VLF 
Ez(l)  - LF 
Ex{l)  - LF 
£*(«)  - LF 
Ez(l ) - BCB 
Bz(l)  • BCB 
BX(Q ) • BCB 
Altimeter 


can  be  replaced  with  SLF 

with  zero 
without  zero 


Not  used 
Not  used 
Not  used 

Flight  number  • remote  entry 

N/C 

Fiducials  lOVFiducials  101/ 

Fiducials  lOVFiducials  lOVFiducials  1C f/ 


Not  used 
Not  used 


• TH  - 10*  hours,  UH  * 1 hour,  TM  - 10*  min,  UM  - 1 min. 
TS  - 10‘  sec,  US  " 1 sec. 


Metrodata.  Each  channel  on  the  tape  recorder  records  dipolar  information  in  a three-digit  format  plus 
sign  information.  Thus,  the  sensitivity  of  the  system  is  one  part  in  999.  All  9 E-phase  channels 
and  the  altimeter  are  recorded  on  the  tape. 

The  tape  recorder  has  a clock  which  is  reset  to  local  real  time  at  the  start  of  a flight  and  incre- 
ments in  seconds.  It  is  internally  formatted  onto  the  tape  and  appears  on  channels  1 and  2.  Chan- 
nel I records  10's  of  hours/unit  hours/10's  of  minutes,  and  channel  2 records  unit  minutes/10*s  of 
seconds/unit  seconds. 

The  fiducial  system  is  based  around  a 5-digit  clock.  Parallel  with  this  clock,  a 5-digit  counter 
is  mounted  m front  of  the  navigator  and  a 5-digit  counter  is  mounted  in  the  camera  to  mark  the  film. 
The  5 digits  of  the  fiducial  information  are  also  recorded  directly  onto  the  tape.  The  entire  system 
is  binary  code  digital  (BCD)  generated.  The  drive  information  for  the  clock  comes  from  an  intervalo- 
meter  which  gives  a clock  pulse  every  1.1  seconds.  Each  pulse  increments  the  clock  by  one  count 
and  cycles  the  tape  recorder.  Every  tenth  pulse  flashes  the  fiducial  count  onto  the  flight  path  cam- 
era film  and  places  a fudicial  on  the  chart  recorder.  Channel  19  is  used  to  indicate  the  flight  num- 
ber, which  is  entered  manually. 

The  navigator  has  the  capability  to  indicate  on  or  off  line  and  to  place  manual  fudicials  on  the 
chart  recorder  and  on  the  flight  path  recovery  film.  The  designations  of  channels  are  listed  in 
Table  II. 
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Figure  17.  Typical  photo  strip  from  flight  camera  on  plane,  from  an  altitude  of  approximately  60  m. 
Landmarks  such  as  rivers  and  roads  were  used  to  place  flight  lines  correctly  on  base  map. 

Flight  path  recovery  caaiera 

A 35-mm  flight  path  camera  manufactured  by  Flight  Research  is  used  to  record  the  flight  path 
of  the  survey  aircraft.  The  camera  is  equipped  with  a fish-eye  lens  and  was  modified  to  give  con- 
tinuous strip  film.  An  example  of  a photo  strip  is  shown  in  Figure  17.  Timing  marks  as  well  as 
the  manual  fiducials  made  by  the  navigator  are  placed  on  the  film. 

Altimeter 

The  height  of  the  aircraft  above  ground  is  determined  by  a Bonzer  radio  altimeter.  The  output 
of  this  altimeter  is  recorded  on  the  analogue  chart  recorder  as  well  as  on  magnetic  tape. 


DATA  REDUCTION 


Horizontal  control 

The  position  of  the  survey  aircraft  along  the  survey  lines  is  determined  from  the  35-mm  flight 
path  recovery  film.  Recognizable  landmarks  on  the  film  are  identified  on  the  photomosaic  or  topo- 
graphic base,  and  the  appropriate  fiducials  are  identified  on  the  analogue  record  and  on  magnetic 
tape.  A straight  flight  path  between  recognized  landmarks  and  a constant  speed  between  fiducial 
points  are  assumed. 

Computation  of  apparent  resistivity 

The  first  step  in  the  calculation  of  the  apparent  resistivity  is  the  determination  of  the  ratio 
£X(Q)/EZ(/)  from  calibration  data  and  flight  data.  During  the  survey  a quick  approximation  of  re- 
sistivity can  often  be  made  from  the  analogue  chart.  In  the  determination  of  EX(Q)/EZ(I) , several 
factors  must  be  taken  into  account: 

1.  The  direction  of  the  flight  line  with  respect  to  the  direction  towards  the  transmitter.  The 
horizontal  component  is  a sinusoidal  function  of  the  coupling  angle  between  the  flight  direction 
and  the  direction  to  the  transmitter.  This  is  illustrated  in  Figure  18,  in  which  data  from  ground- 
based  measurements  of  in-phase  Ex(/)  and  quadrature-phase  EX(Q)  are  plotted  as  a function  of 
compass  bearing  (instrument  magnetic  bearing).  Since  the  horizontal  antenna  is  at  90°  to  the  flight 
direction,  maximum  coupling  occurs  when  the  angle  between  flight  line  and  the  direction  to  the 
transmitter  is  90°. 
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0 40  80  120  160  ZOO  240  260  320  360 


Instrument  Magnetic  Bearing,  degrees 

Figure  18.  Ground  measurements  -of  the  in-phase  and  quadrature-phase  com- 
ponents of  W as  a function  of  compass  bearing- 

2.  The  gain  settings  during  calibration  and  operation. 

3.  The  field  strength  of  the  vertical  electric  field  Ey  during  calibration  and  operation. 

Taking  all  these  factors  into  account,  the  ratio  Ex(QVEy  is  computed  from  the  following 
relation: 

E*«?>  0.34  A cal  , Ezcal  E*W° Per 

x x * ltd) 

E z Ex(Q)cal  A oper  Eloper  sinn 

where  Ex(Q)cal  -=  the  quadrature-phase  component  of  Ex  at  calibration  (in  millimeters  on  the  analogue 
record) 

A cal  - attenuator  setting  during  calibration 
i4  opet  = attenuator  setting  during  survey 

Ez  cal  “ Ey  at  time  of  calibration  of  EX(Q)  cal  (in  millimeters  on  the  analogue  record) 

Ez  oper  Ez  during  the  survey  at  the  position  where  the  resistivity  is  being  calculated 
(in  millimeters  on  the  analogue  record) 

EX(Q)  oper  EX(Q)  at  the  position  where  the  resistivity  is  calculated  (in  millimeters  on  the 
analogue  record) 

a the  angle  between  the  average  true  azimuth  of  a flight  line  and  the  true  azimuth 
from  the  center  of  the  survey  area  to  the  transmitter.  In  case  of  maximum  coupling 
90°  (sin  90°  1). 

These  quantities  are  registered  as  volts  on  the  magnetic  tape  or,  when  determined  from  the  analogue 
records,  in  millimeters.  When  the  transmitters  are  close  to  the  survey  area,  as  was  the  case  m some 
BCB  and  LF  surveys,  the  value  of  sin  a may  vary  along  the  flight  line.  This  is  taken  into  consid- 
eration during  computation;  the  true  azimuth  of  the  flight  line  is  determined  between  fiducial  points 
and  a new  « is  determined  for  each  section  of  the  flight  line,  but  between  fiducials  a is  assumed 
constant.  For  a given  calibration  and  a,  eq  10  reduces  to 

EX(Q)  E (Q)  oper 

_1_  _ c -J (15) 

Ez  Ez  oper 


where 


Table  M.  Survey  specifications  aad  transmitter  data. 
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c -Off—  x <r!1.  x ELc!!.  . (16) 

EX(Q ) cal  * A oper  sin  a 

The  apparent  resistivity  is  subsequently  computed  from  eq  11. 

In  Table  III  the  specifications  of  the  various  study  areas  discussed  in  this  report  are  listed. 

Computer  processing  of  data 

The  computed  apparent  resistivities  at  each  fiducial  point  are  listed  in  the  form  of  resistivity 
profiles  for  each  survey  line  (P-plot).  The  fiducial  points  are  given  on  the  vertical  axis  and  the 
three  resistivities  are  printed  above  the  corresponding  fiducial  points  using  three  different  symbols. 
The  vertical  scale  of  the  fiducial  marks  is  arbitrary.  Whenever  the  calculated  value  of  pa  is  nega- 
tive - which  may  be  due  to  powerline  interference,  scattering,  and  rero  error  - negative  values 
are  plotted  on  the  printouts. 

Plotting  and  contouring  of  data 

The  fiducial  points  common  to  the  base  maps  and  magnetic  tape  are  merged  and  the  computed 
resistivity  data  are  plotted  along  the  flight  lines,  prom  these  plots  the  data  are  contoured  by  hand. 
Examples  of  resistivity  profiles  along  a flight  line  are  given  in  Figure  19. 

Filtering  of  B-phase  data 

To  smooth  the  observed  vertical  electric  field,  a low  pass  filter  is  applied  to  the  field  before 
computing  the  apparent  resistivity;  the  computed  apparent  resistivities  are  also  filtered.  The  ver- 
tical electric  field  data  and  the  resistivity  data  were  filtered  for  each  survey  area. 
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PROBLEM  AREAS  OP  THE  E-PHASE  TECHNIQUE 

The  three  major  problems  with  the  system  at  present  are  horizontal  control  (navigation),  zero 
errors  and  interference. 

Horizontal  control 

The  flight  lines  are  located  on  the  map  by  merging  recognizable  landmarks  on  the  film  strip 
with  flight  time  locations  marked  by  the  navigator  on  the  topographic  or  mosaic  base  at  the  time  of 
the  survey.  The  error  in  this  technique  clearly  depends  much  on  the  character  of  the  terrain:  when 
a few  landmarks  stand  out  as,  for  example,  in  the  Moose  Creek  Dam  area,  the  navigator  has  great 
difficulty  in  spacing  the  flight  lines.  In  the  Chena  Dam  area  flight  lines  often  cross.  The  error  in 
the  horizontal  control,  therefore,  varies  much  from  one  area  to  another.  In  these  surveys  no  effort 
was  made  to  set  out  special  markers  or  navigational  beacons  to  help  reduce  the  error  in  horizontal 
control.  This  problem  naturally  increases  greatly  with  very  closely  spaced  flight  lines.  The  sur- 
veyed areas  were  as  a general  rule  flown  at  Vlt*mile  (0.16-km)  intervals. 

Zero  error 

Calibration  is  accomplished  at  the  beginning  and  end  of  a survey  by  tilting  the  antennas,  and 
this  calibration  is  assumed  constant  during  the  survey.  Drift  in  the  instrument  is  very  small,  but 
Other  factors  affecting  calibration,  such  as  the  direction  of  £z(f)  and  EX(Q),  may  change  due  to 
scattering  of  the  field  and  mixed-path  propagation.  Negative  resistivities  did  occur  on  approxi- 
mately 1%  of  the  total  line  mileage  and  they  were  probably  partly  due  to  zero  error.  One  way  to 
reduce  zero  errors  is  to  use  a continuously  rotating  antenna,  so  that  a zero  is  obtained  and  the 
wave  tilt  is  calibrated  at  each  point  in  the  survey. 

Iztorfermee 

Interference  is  encountered  1)  whenever  the  flight  path  encounters  man-made  objects  that  alter 
the  local  electric  field,  and  2)  whenever  topographic  features  diffract  or  scatter  the  electric  field. 
The  interference  caused  by  man-made  objects  is  easily  recognized  and  results  in  large  anomalies 
in  the  records.  The  data  in  such  instances  are  discarded.  The  interference  caused  by  topographic 
features  is  more  difficult  to  recognize  and  the  effects  of  these  features  must  await  further  studies. 
In  general,  no  systematic  variation  at  the  apparent  resistivity  with  topography  was  noticed,  but 
further  investigations  on  this  aspect  are  planned. 


GROUND  CONTROL  IN  STUDY  AREAS 

The  sites  were  selected  because  of  the  availability  of  ground  control.  The  amount  of  ground 
control  available  consisted  of  geological  maps  prepared  by  the  U.S.  Geological  Survey  (USGS)  and 
soils  maps  published  by  the  U.S.  Department  of  Agriculture;  in  addition,  in  some  areas  the  Alaskan 
District,  Corps  of  Engineers,  provided  results  from  the  preliminary  site  surveys  for  the  Chena  Flood 
Control  Project. 

The  control  information  varied  from  one  site  to  another.  In  some  cases,  no  detailed  surficial 
geological  mapping  was  available  and  control  was  based  on  information  from  widely  spaced  drill 
logs.  The  study  areas  and  the  type  of  control  available  in  them  are  listed  in  Table  IV. 

The  geological  unitt  used  in  developing  the  comparisons  between  surficial  geology  and  re- 
sistivity were  directly  taken  from  the  geological  maps  prepared  by  the  USGS;  the  characteristics 
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of  the  units  are  listed  in  Table  V (Pdwd  1958.  Williams  et  al.  1959).  In  order  to  use  the  soil 
map  (Rieger  et  al.  1963).  soil  units  were  grouped  in  larger  classes;  this  grouping  will  be  discussed 
later. 


Table  IV.  Sanmary  of  ground  control  available  in  each  survey  area. 

Study  area  Control 


Goldstream  area 
Site  2 

Chena  Hot  Springs  Road 
Moose  Creek  Dam  area 


Surficial  geology  maps,  soils  maps,  and  drill  logs 
Surficial  geology  maps,  soils  maps,  and  drill  logs 
Drill  logs 

Geology  and  limited  drill  logs 


Table  V.  Geological  units  from  Fairbanks,  Alaska,  D-l  and  D-2  Geological  Quadrangles. 


Unit 

Material  type 

Thickness 

Permafrost 

Qsu 

Perennially  frozen 

undifferentiated 

silt 

Silt-organic 
ground  ice 

1 to  30  m (D-i) 

Depth  to  permafrost  0.45  to 
1.20  n>  on  lower  slopes 
and  creek  valley  bottoms, 

1.5  to  6.10  m near  contacts 
with  Of  and  pCbc.  Active 
layer  0.45  to  i,20  m.  Perma- 
frost 1 to  52  m. 

Qso 

Perennially  frozen 
organic  silt 

Organic-silt 
ground  ice 

15  to  > 30  m (D-2) 

Active  layer  0.45  to  1 m. 
Permafrost  at  least  30 
m thick. 

Qs 

Flood  plain  swale 
and  slough 
deposits 

Silt  and  organic 
material 

Generally  < 3 m 
Rarely  3 to  7.6  m 

Permafrost  depth  0.45  to 
6.1m.  Permafrost  1.5  to 
9.1  m thick. 

Of 

Fairbanks 

silt 

Low  density 
silt 

Mapped  when  • 1 m 
on  hill  tops,  and 
as  much  as  30  in 
on  midslopes 

No  permafrost . 

Qg 

Creek  gravel 
incl  tailings 

Gravel  with 
local  fine  lenses 

1 to  76  m 

Only  locally  frozen. 

Qal 

Flood  plain 
alluvium 

Silts  over  sand 
and  gravel 

Silt  0.3  to  4,6  m 
thick  over  sand 
and  gravel  of 
varying  thickness 
up  to  approx  213  m 
near  the  Tanana  River. 

Permafrost  discontinuous, 
unfrozen  lenses,  layers 
and  vertical  zones. 
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Unit 

pGbc 

Bedrock  Birch 
Creek  schist 


t?r 

Granite 


Table  V.  Continued 

Material  type  Thickness 


Less  than  1 m of 
silt  over  rock 


Less  than  1 tn 
of  silt  over  rock 


Permafrost 


Permafrost  generally  not 
found  except  under  Qsu 
and  Qg  in  valley  bottoms. 


COMPUTER  MODELING  OF  RESISTIVITY  PROFILES  IN  CENTRAL  ALASKA 

Natural  ground  usually  consists  of  several  layers  of  different  resistivities  and  the  area  in 
Fairbanks,  Alaska,  is  no  exception.  Typical  resistivity  profiles  of  the  major  surficial  geology 
units  found  in  the  vicinity  of  Fairbanks,  Alaska,  are  schematically  listed  in  Table  VI.  Seasonal 
changes  in  the  resistivity  profile  occur  in  the  top  3 m of  soil  and  are  most  severe  when  an  organic 
surface  layer  is  present.;  The  variations  in  thickness  of  the  units  listed  in  Table  VI  correspond  to 
the  variations  suggested  for  the  surficial  geological  units  by  the  local  geological  studies.  The 
ranges  of  resistivity  listed  are  based  on  laboratory  and  field  measurements.  The  results  of  the 
laboratory  measurements  were  given  in  section  Resistivities  of  earth  materials  (p.  3).  Some  of  the 
field  measurements  are  discussed  next. 

Galvanic  resistivity  measurements  (Keller  and  Frischknecht  1966)  were  made  with  the  elec- 
trodes in  a Wenner  array  on  each  of  the  units  listed  in  Table  V.  In  the  Wenner  method,  the  spac- 
ing between  the  electrodes  {"a  spacing")  is  plotted  versus  the  apparent  resistivity. 

Increasing  the  spacing  between  the  electrodes  increases  the  depth  to  which  subsurface  infor- 
mation may  be  obtained.  The  depth  of  penetration  is  about  equal  to  the  distance  between  the 
electrodes.  When  the  earth  is  not  homogeneous,  but  has  several  layers  of  different  resistivities, 
Pj.’Pg,  p3...,  and  of  thicknesses,  dj,  dg,  dg..  . the  value  of  the  apparent  resistivity  changes  with 
electrode  spacing.  In  the  field  ptt  is  measured,  and  pj,  pg,  Pg,...  and  dj,  dg,  d3>...  are  obtained 
by  fitting  a computer  solution  to  the  field  data. 

In  Figures  20a-d  the  solid  lines  represent  the  computed  curves  for  the  ground  conditions 
indicated  on  the  figures.  In  the  computer  program,  the  apparent  resistivity  obtained  in  a Wenner 
electrode  configuration  can  be  simulated  for  ground  of  up  to  four  layers.  A best  fit  is  obtained  by 
successive  trial  inputs  of  layer  resistivities  and  thicknesses. 

Figures  goa  and  b are  representative  of  Wenner  spreads  on  the  units  Qsu  and  Qso  under 
summer  conditions.  In  Figure  20&  the  Wenner  spread  is  shown  for  a thawed  organic  mat,  0.3  m 
thick  and  with  a resistivity  of  140  ohm-m,  overlying  permafrost  with  a resistivity  of  7000  ohm-m. 

In  this  case  the  spacing  between  the  electrodes  was  not  opened  far  enough  to  indicate  the  bottom 
of  the  permafrost.  Because  of  the  high  resistivity  of  the  permafrost  layer,  the  frozen  ground  most 
likely  has  a high  ice  content.  In  Figure  20b  two  Wenner  spreads  are  shown  for  permafrost  sites 
where  the  resistivity  of  the  frozen  ground  is  between  500  and  700  ohm-m  and  the  depth  to  the  per- 
mafrost table  is  deeper. 
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Table  VL  Layering  and  resistivity  range*  expected  on  ■ajar  ssrilclal  geological  utte. 


Thickness 

Resistivity  (ohm-m) 

[feit 

Layer 

(■) 

Water  Sommer 

Qsu,  Qso 

Active  layer 

0.2  to  3 

> 1,000  60  to  1&0 

Unfrozen  silt 

0 to  6 

100  to  800 

Frozen  silt 

5 to  60 

600  to  10,000 

Bedrock  (schist) 

Basement* 

100  to  > 3,000 

or 

Seasonal  frost  layer 

t to  3 

> 1,000  100  to  800 

Unfrozen  silt 

t to  30 

100  to  800 

Bedrock  (schist) 

Basement 

100  to  2,000 

Qal 

Seasonal  frost  layer 

1 to  3 

> 1,000  100  to  1,000 

Unfrozen  silt 

0 to  7 

100  to  800 

Gravel 

0 to  800 

1,000  to  10,000 

Bedrock  (schist) 

Basement 

100  to  2,000 

• Lower  layer  of  sufficient  depth  so  that  it  can  be  considered  infinite  at  VLF  radiation. 


Typical  Wenner  spreads  on  the  units  Of  are  shown  in  Figure  20c.  These  units  are  most  fre- 
quently situated  on  the  south-facing  hillsides  and  have  in  general  a silt  cover  less  than  30  m 
thick;  the  apparent  resistivities  of  these  units  at  “a"  spacings  (separation  distance  between 
electrodes  in  a Wenner  array)  larger  than  30  m are,  therefore,  mainly  determined  by  bedrock 
resistivity. 

Figures  20d  ana  e are  typical  examples  of  Wenner  spreads  on  the  Qal  units.  Curve  A of  Fig- 
ure 20d  represents  virtually  no  silt  cover  over  the  gravel;  curve  B represents  an  unfrozen  silt  cover 
about  2.5  m thick.  An  example  of  a Wenner  spread  on  a unit  Qal  in  the  winter  is  shown  in  Figure 
20e.  The  seasonal  frost  is  evidenced  by  a layer  of  high  resistivity  dj  at  the  surface,  followed  by 
a layer  of  unfrozen  silt  and  gravel  dg. 

Table  VI  and  the  Wenner  spreads  demonstrate  the  large  variation  in  resistivity  that  can  occur 
in  the  surricial  geological  units  of  the  D*1  and  D-2  Fairbanks,  Alaska,  Geological  Quadrangles. 

The  ranges  in  values  are  due  both  to  variation  in  material  type  (e.g.  ice  content)  and  to  variation 
in  the  thickness  of  each  layer. 

Using  the  limits  of  the  resistivities  and  thicknesses  of  each  unit  listed  in  Table  VI,  the  vari- 
ation in  the  apparent  resistivity  expected  from  an  E-phase  survey  can  be  computed  from  eq  11  at 
each  desired  frequency. 

In  the  units  Qsu  and  Qso,  the  main  variations  at  VLF  are  caused  by  the  thickness  of  the  silt 
cover,  the  resistivity  of  the  basement,  and  the  thickness  and  resistivity  of  the  active  layer.  In 
Figure  21  the  computed  apparent  resistivity  at  VLF  and  LF  in  the  units  Qsu  and  Qso  is  given  as 
a function  of  the  thickness  of  frozen  silt  cover  d j at  two  values  of  the  resistivity  of  the  basement 
(p4  = 1000  ohm-m  and  p4  «=  300  ohm-m)  in  the  presence  of  an  active  layer.  The  apparent  resistivity 
at  VLF  in  the  unit  Qsu  is  bounded  by  a low  of  400  ohm-m  and  a high  in  excess  of  2500  ohm-m.  The 
variation  in  apparent  resistivity  is  lesR  at  LF  because  of  the  screening  eff  jet  of  the  conductive 
surface  layer.  Also,  because  the  thickness  of  the  silt  cover  in  the  unit  Qso  is  often  in  excess  of 
15  m,  a narrower  resistivity  range  is  expected  for  Qso  than  for  Qsu. 

The  expected  resistivity  variation  in  the  unit  Qf  is  more  difficult  to  model  because  of  the 
limited  amount  of  silt  cover  over  bedrock.  Because  of  the  thin  (<  5-m)  silt  covet  over  bedrock, 
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Figure  21.  The  computed  apparent  resistivi- 
ties for  the  E-phase  method  at  VLF  and  LF  on 
the  units  Qso  anJ  Qsu  as  a tunction  of  thick- 
ness of  the  frozen  silt  cover. 


Figure  22.  The  computed  apparent  resistivity 
at  VLF  typical  for  the  unit  Qal.  The  thick- 
ness of  the  silt  cover  dj  is  varied  from  5 to  10 
m and  pa  is  plotted  as  a function  of  the  thick- 
ness of  gravel  dg.  The  bedrock  ps  is  assumed 
conductive  at  50  ohm-m. 


Figure  23.  The  computed  apparent  resistivity 
at  LF  typical  for  the  unit  Qal.  The  resistiv- 
ity is  plotted  as  a (unction  of  the  thickness  of 
the  silt  ewer  dj. 


Table  VII.  VLF  (NLK,  18.6  kHz)  ground-baaed  apparent  resistivity 
measurements  os  roek  sites  in  the  Ynkon-Tanana  Uplands. 


Site 

Apparent 

resistivity 

(ohm-m) 

Remarks 

1 

160  to  180 

Schist,  severely  decomposed 

2 

130  to  160 

Schist,  severely  decomposed 

3 

400  to  500 

Schist,  badly  weathered 

4 

600  to  1100 

Schist,  badly  fractured 

& 

2000 

Intrusive  rock,  weathered  and  fractured 

e 

4000  to  5000 

Oranitic  intrusive  rock,  jointed  and  fractured 

the  VLF  values  are  mainly  de*ermined  by  bedrock  resistivities.  Table  VII  shows  VLF  ground-based 
measurements  on  rock  outcrops  in  the  Yukon  Tanana  Uplands  (Sellmann  et  al.  1974).  Since  some 
c*  these  values  ate  from  rock  quarries  where  the  quality  of  rock  may  perhaps  be  better  than  nor- 
mally encountered,  VLF  apparent  resistivities  on  Of  probably  fall  between  50  and  1000  ohm-m. 
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The  LF  resistivity  values  should  correspond  to  the  apparent  resistivities  obtained  with  the 
Wenner  method  at  spacings  between  10  and  20  m.  In  the  summer  these  values  vary  from  500  ohm*m 
to  2000  ohm-m  (see  Fig.  20c),  and  may  be  higher  in  the  winter. 

In  the  Qal  unit,  the  thickness  of  the  unfrozen  silt  cover,  the  thickness  and  resistivity  of  the 
flood  plain  gravels,  and  the  resistivity  of  the  basement  are  the  main  parameters  influencing  pA  at 
VLF.  Figure  22  shows  the  dependence  of  the  computed  apparent  resistivity  at  VLF  on  two  of 
these  factors:  the  thickness  of  the  siU  cover  and  the  thickness  of  gravel.  In  Figure  23  the  com- 
puted apparent  resistivity  ii  given  at  LF  as  a function  of  the  thickness  of  the  silt  cover  over 
gravel.  The  resistivity  of  the  basement  has  no  effect  in  the  situation  shown  because  of  the  thick 
intermediate  gravel  layer. 

The  major  findings  from  the  computer  modeling  studies  as  they  apply  to  delineating  the  surfi- 
cial  geological  units  of  the  Fairbanks,  Alaska,  D*1  and  D-2  Geological  Quadrangles  are  summarized 
below: 

1.  Apparent  resistivities  between  400  ohm-m  and  > 2500  ohm-m  are  expected  at  VLF  on  the 
units  Qso  and  Qsu.  The  resistivities  increase  with  the  depth  of  permafrost.  The  VLF  resistivi- 
ties on  the  units  Qf  are  expected  to  be  less  than  500  ohm-m;  this  expectation  is  mainly  based  on 
ground  measurements. 

2.  In  the  units  Qso  and  Qsu.  the  apparent  resistivity  at  VLF  is  greater  than  the  value  at  LF 
because  of  the  screening  effect  of  the  conductive  surface  layer.  Based  on  data  from  Wenner 
spreads,  the  values  at  LF  are  greater  than  700  ohm-m  on  the  unit  Qf. 

3.  In  the  unit  Qal,  apparent  resistivities  are  mainly  determined  by  the  thickness  of  the  silt 
cover. 


RESULTS 

Although  the  resistivity  data  were  often  contoured  at  50  ohm-m,  all  maps  in  this  report  have 
only  3 contour  lines:  250,  500  and  1000  ohm-m.  This  provides  four  resistivity  ranges,  which  were 
found  to  correspond  well  to  the  natural  limits  expected  for  various  ground  conditions,  permafrost 
and  geological  settings. 

A problem  in  comparing  geological  maps  with  apparent  resistivity  data  is  the  large  degree  of 
variation  in  resistivity  allowed  within  a geological  unit,  as  illustrated  in  Figures  21  and  22.  Fro- 
zen silt  covers  in  the  units  Qso  and  Qsu  range  in  thickness  (tom  1 m to  60  m.  Even  though  these 
variations  exist,  the  overall  differences  often  appear  significant  enough  to  allow  a distinction  be- 
tween units  on  the  basis  of  resistivity.  This  type  of  approach  to  evaluating  the  data  was  also 
considered  in  keeping  with  the  intended  application  of  the  airborne  data  as  a regional  tool. 

Goldstreaa  site 

The  Goldstreaui  area  has  a geological  setting  typical  of  many  stream  valleys  in  the  uplands 
between  Fairbanks  and  the  Yukon  Valley;  in  it  are  found  several  common  geologic  and  permafrost 
associations: 

1.  Silt  units,  which  are  further  subdivided  into: 

(a)  thawed  silt  of  wind-blown  origin  common  on  south-facing  exposures,  indicated  by  Qf 

(b)  frozen  silt  units  common  to  north-facing  hills,  indicated  by  Qsu 

(c)  frozen  silt  units  common  to  valley  bottoms,  indicated  by  Qsu  or  Qso 

2.  Gravel  dnits,  which  are  associated  with: 

(a)  contemporary  streams 

(b)  old  channels,  which  may  be  found  reworked  in  local  tailing  piles.  The  older  gravel 
deposits  may  be  perennially  frozen.  Gravel  units  are  indicated  by  Qg. 


' the  Fairbanks,  Alaska, 


Elevation,  ft 
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3.  Bedrock  units.  Bedrock  is  generally  covered  by  silt  or  gravel,  but  may  be  exposed  on  steep 
slopes  and  on  some  ridge  tops.  The  common  bedrock  is  metamorphic,  often  weathered,  and 
badly  decomposed;  but  intrusive  rock  occurs  locally.  Exposed  bedrock  sites  with  Birch 
Creek  schist  are  indicated  by  pebc  and  when  the  exposed  rock  is  intrusive  by  gr. 

The  VLF  and  LF  resistivity  contour  maps  with  the  surficial  geological  map  of  the  D-l  Fair* 
banks  Quadrangle  (Pdwd  1958)  are  given  in  Figures  24  and  25.  respectively. 


VLF  resistivities.  In  comparing  the  contoured  VLF  data  with  the  geological  map,  a noticeable 
degree  of  correlation  is  evident.  Significant  differences  are  observed  between  the  resistivities 
of  the  thawed  units  Of  and  those  of  the  perennially  frozen  silt  units  Qsu  and  Qso. 

The  western  half  of  the  Goldstream  study  area  is  largely  a broad  valley  with  an  extensive  silt 
cover.  The  silts  extend  across  the  valley  from  ridge  top  to  ridge  top.  The  perennially  frozen  silts 
are  common  in  the  valley  bottoms  and  north*facing  hillsides,  and  grade  into  the  thawed  silts  on  the 
more  south-facing  exposures.  Figure  23  is  a typical  cross  section  of  the  valley  along  A*A';  it  il* 
lustrates  the  general  relation  between  the  geological  units.  The  cross  section  is  based  on  the 
topographic  map  and  the  results  of  7 drill  holes  (Pdwd  1958). 

The  frozen  silts  in  the  valley  (Qso,  Qsu)  are  characterized  by  high  resistivities,  often  in  ex- 
cess of  1000  ohm-m.  The  values  decrease  at  the  valley  margins,  where  the  thawed  silt  units  com- 
monly have  apparent  resistivities  of  less  than  500  ohm-m,  and  often  below  250  ohm-m.  The  500-ohm-m 
contour  is  often  near  the  mapped  surficial,  geological  boundary  of  the  units  Qf  and  Qsu. 

The  QT  units  at  the  northern  margin  of  the  valley  are  interrupted  by  a unit  of  Qsu  situated  in  a 
small  tributary  of  Goldstream.  The  resistivity  values  also  reflect  this  interruption  with  a break  from 
the  low  resistivities  to  the  higher  values  indicative  of  perennially  frozen  ground  in  the  tributary 
valley.  Several  other  smaller  low  resistivity  zones  are  found  associated  with  the  Qf  unit  on  the 
southern  valley  margin. 


Resistivities  in  excess  of  500  ohm-m  are  not  restricted  to  frozen  silt  sections,  but  are  also 
found  associated  \M«h  gravel  hi  the  northern  part  of  the  study  area  and  often  with  bedrock  on  the 
northwestern  valley  margin. 
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Figure  26.  A typical  cross  section  of  Goldstream  Valley  along  A- A'  shown  in  Figure  24. 
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There  are  also  zones  of  apparent  disagreement,  such  as  the  small  zones  of  low  resistivities  in 
the  valley  bottoms  mapped  as  Qsu.  In  some  cases,  it  is  possible  that  the  resistivity  data  are  more 
accurate  than  the  contour  map  in  terms  of  permafrost,  particularly  since  permafrost  is  subject  to 
change  in  areas  of  human  activity  such  as  the  Fairbanks  area.  Only  more  detailed  work  will  permit 
explanations  of  these  anomalies  found  throughout  the  study  area. 

In  summary,  the  VLF  map  in  Goldstream  Valley  shows  in  general  a good  correlation  between 
the  apparent  resistivities  and  the  mapped  geological  boundaries  of  the  frozen  and  thawed  silt  units. 
The  major  problem  in  making  these  ground  correlations  is  the  extreme  variation  within  the  geological 
units  discussed  in  the  preceding  section. 

LF  resistivities.  When  the  LF  map  is  compared  with  the  VLF  map  several  differences  are 
observed: 

1.  In  the  unit  Of  the  resistivities  at  LF  are  higher  than  those  at  VLF;  whereas  at  VLF  values 
are  in  general  lower  than  500  ohm-m,  at  LF  in  several  sections  resistivities  are  higher  than  1000 
ohm-m. 

2.  In  the  perennially  frozen  units  Qsu  and  Qso,  the  trend  is  reversed:  the  resistivities  at  LF 
are  lower  than  those  at  VLF. 


To  investigate  these  resistivity  trends  at  VLF  and  LF  further,  normalized  distribution  histo- 
grams were  made  by  tabulating  every  other  resistivity  value  of  a P-plot  with  the  geological  units 
under  the  flight  path.  The  data  were  stored  in  computer  files  and  the  resistivity  distribution  curves 
were  made  from  the  files  with  a sorting  program.  Figure  27  gives  the  normalized  distribution  histo- 
grams for  the  units  Qf,  Qso  and  Qsu  at  LF  and  VLF  resistivities,  respectively.  In  addition  to  the 
resistivity  intervals  0-250,  250-500  and  500-1000  ohm-m  used  on  the  maps,  the  intervals  1000-2000 
and  greater  than  2000  ohm-m  are  given.  These  distribution  curves  confirm  the  remarks  made  in  the 
preceding  discussion:  1)  at  VLF  the  units  Qso  and  Qsu  have  a higher  density  of  resistivities  in 
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Figure  27.  Normalized  distribution  histograms  of  resistivities  on  the 
units  of  Qf,  Qsu  and  Qso. 
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Table  Vm.  bleaa  apparent  resistivities  ia  several  geological  —its. 


Goldsuetm  area 


Geological 

p.i  ohm-m 

unit 

VLF,  18.6  kl lx 

LF.  876  kS * 

Qso,  high  i ce 

1538 

702 

Qsu 

1183 

960 

Qf 

897 

1145 

Site  2 

pt,  ohm-a  “ 

VLF,  18.6  kHz LF,  S57  kHz 


711  300 

345  424 


excess  of  1000  ohm-m  than  the  unit  Of.  2)  At  LP  the  density  of  resistivities  in  excess  of  1000  ohm*m 
decreases  for  the  units  Qsu  and  Qso  and  increases  for  the  unit  Of  compared  with  the  density  of  re* 
sistivities  at  VLF.  This  reversal  is  further  confirmed  by  the  mean  values  of  resistivity  listed  in 
Table  VIII. 

In  Figure  28  the  apparent  resistivities  of  the  units  Qsu,  Qso  and  Of  are  plotted  as  a function  of 
frequency.  A value  for  BOB  was  obtained  from  sites  where  the  quality  of  the  data  was  good.  The 
shaded  area  indicates  the  range  of  values  observed  and  the  mean  values  fall  within  this  range.  The 
dotted  line  gives  the  computed  values  for  the  ground  conditions  indicated.  Thus,  in  the  frozen  units 
Qsu  and  Qso,  where  pa  decreases  with  frequency,  a conductive  surface  layer  d j and  d2  overlies  a 
layer  d3  of  high  resistivity  p3,  followed  by  a basement  resistivity  p4  of  300  ohm-m.  The  computed 
value  of  pa  is  sensitive  to  all  the  parameters  in  the  model.  For  example,  if  the  thickness  of  frozen 
silt  is  increased,  the  value  of  pa  at  VLF  becomes  larger.  Also,  the  resistivity  p3  must  be  in  ex- 
cess of  5000  ohm-m  to  yield  the  strong  decrease  in  resistivity  with  frequency  observed  in  the  data. 

Since  the  survey  was  made  in  early  spring,  when  there  was  only  a very  Bmall  amount  of  thaw, 
the  presence  of  a conductive  layer  at  the  surface  was  initially  difficult  to  explain.  Later  laboratory 
work  on  samples  from  the  organic  surface  layer  common  to  permafrost  showed  that  the  resistivities 
decrease  rapidly  when  the  temperature  approaches  0°C. 
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Figure  28.  The  apparent  resistivity  measured  as  a function  of  frequency  on  the  geological  units 
Qsu,  Qso  and  Qf,  The  shaded  area  gives  the  range  of  values  encountered.  The  dotted  line  gives 
the  computed  values  for  the  ground  condition  indicated. 
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In  summary,  the  multifrequency  airborne  resistivity  survey  in  the  Goldstream  area  has  shown 
that  the  VLF  resistivity  map  shows  good  correlation  between  the  mapped  boundaries  of  frozen  and 
thawed  silt  units.  The  correlation  between  resistivity  at  LF  and  surficial  geology  is  less  evident. 
Large  resistivity  changes  occur  between  VLF  and  LF;  on  frozen  units  the  resistivity  decreases 
with  increasing  frequency,  and  on  thawed  units  the  opposite  is  the  case. 

Site  2 

This  study  area  (Fig.  1)  is  situated  on  the  northeast  of  Fairbanks  and  has  a general  N*S  trend. 

It  covers  a swath  that  makes  the  transition  from  the  uplands  including  Gilmore  Dome  into  the  low 
flood  plains  of  the  Tanana  Valley.  The  surficial  geology  includes  the  fine-grained  soils  of  the  up- 
lands, some  of  which  are  perennially  frozen,  and  the  mixed  alluvial  materia!  in  the  Tanana  Valley. 
Also,  bedrock  types  are  variable,  including  both  metamorphic  and  igneous  rock  types. 

The  area  was  flown  4t  LF,  VLF  and  BOB;  both  detailed  surficial  geological  and  soil  maps  are 
available  for  this  area,  providing  considerable  detailed  control  for  assessing  the  results. 

The  VLF,  LF  and  BCB  resistivity  maps  of  site  2 are  shown  in  Figure  29  with  a geological  map 
and  in  Figure  30  with  a soils  map. 

VLF  resistivities.  The  VLF  data  again  provide  the  most  obvious  correlation  between  the  per- 
ennially frozen  units  Qsu  and  Qso,  and  the  thawed  unit  Of.  The  central  part  of  the  study  area  be- 
tween the  Chena  River  and  the  uplands  is  mapped  as  Qsu.  This  area  also  has  consistently  high 
resistivity  values.  The  500-ohm-m  contour  has  the  same  general  limits  as  the  mapped  Qsu.  Im- 
mediately to  the  northeast  is  a mapped  unit  of  Of,  which  has  consistently  low  values  of  resistivity. 
However,  the  adjacent  Qsu  unit  to  the  north  contains  no  resistivity  highs  as  would  be  expected. 

This  zone  and  others  like  it  will  be  subjected  to  further  fV* J work  to  determine  the  actual  condi- 
tions. There  is  a possibility  that  the  mapped  limits  of  some  of  these  units  have  changed  or  are  in 
error.  The  resistivity  distribution  histograms  for  tie  units  Of,  Qsu  and  Qso  at  VLF  and  LF 
shown  in  Figures  31a  and  b,  respectively,  again  substantiate  the  general  conclusion  that  the  re- 
sistivities are  higher  in  units  Qsu  than  in  units  Of,  and  the  mean  values  listed  in  Table  VIII  also 
reflect  this  trend. 

In  the  southern  end  of  this  study  area,  the  surficial  geological  units  associated  with  the  Tanana 
flood  plain  are  found.  There  is  no  apparent  correlation  between  resistivity  contours  and  the  mapped 
boundaries  of  the  units  Qal  and  Qs  and  perhaps  little  ought  to  be  expected.  The  units  Qal  consist 
of  a silt  layer  1 to  5 m thick  over  gravel,  and  the  units  Qs  are  in  old  river  channels,  with  probably 
similar  material  types.  What  the  resistivity  map  might  indicate  is  the  thickness  of  silt  over  gravel, 
but  the  surficial  geology  map  provides  no  control  on  that. 

LF  and  BCB  resistivities.  When  the  LF  and  BCB  resistivity  data  are  compared  with  VLF  data, 
the  same  reversal  in  resistivity  with  frequency  observed  in  the  Goldstream  area  is  found:  in  the 
unit  Qsu  the  resistivity  decreases  with  frequency,  and  in  the  unit  Of  the  resistivity  increases  with 
frequency.  The  resistivity  distribution  histograms  shown  at  LF  in  Figure  31b  and  the  mean  values 
of  resistivity  listed  in  Table  VIII  again  confirm  this  conclusion.  However,  the  reversal  in  this 
area  is  not  as  pronounced  as  in  the  Goldstream  area.  There  are  no  apparent  correlations  between 
the  resistivity  contours  and  the  surficial  geology  at  LF  and  BCB. 

Some  trends  between  resistivity  at  LF  and  BCB  and  data  from  the  Fairbanks  soils  map  can  be 
discerned  from  the  resistivity  contour  map  w:t’.i  the  soil  map  in  Figure  30.  The  soils  were  grouped 
in  classes  according  to  permafrost  distribution  and  thickness  of  active  layer.  Table  IX  lists  the 
classifications  used  and  their  characteristics.  Classes  Ab,  Ta  and  Tb  include  nonpermafrost  soils, 
and  soils  with  permafrost  at  a depth  greater  than  4.6  m.  These  soils  are  associated  with  the  higher 
resistivities  in  the  central  and  northern  parts  of  the  area.  The  soil  units  with  the  permafrost  table 
within  about  1 m from  the  surface,  Aa  and  Pa,  have  resistivities  less  than  500  ohm-m  at  LF  and 
less  than  250  ohm-m  at  BCB. 


UQL 


apparent  resistivity  contour  map  of  site  2,  with  soils  map.  Blue  = <250,  white  = 250-<500,  yel- 
low = 5 00-1000,  and  red  >1000  ohm-m. 
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Figure  31.  Normalized  distribution  histograms  of  resistivities  on 
the  units  Ql,  Qsu  and  Qso. 


Table  IX.  Qrosplag  of  the  soil  salts  from  the  Fairhaaks,  Alaska,  soil  suvey  (Rieger  et  aL  19€3). 


Data 
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analysis 

soil  survey 

classification 

classification 
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Ta 

Ab 

Sc 
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Pa 
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Su 

LP 
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Pb 

Mn 

Ta 

Or 

Qm 

Fm 

Tb 

Fa 

Me 

Me 

Soil  name 

Bradway  fine  sandy  loam 
Tanana  silt  loam 
Salchaket  fine  sandy  loam 
Salchaket  fine  sandy  loam 
Coldstream  silt  loam 
Saulich  silt  loam 
Lemetapeat 
Ester  silt  loam 
Minto  silt  loam 
Qilmore  silt  loam 
Qilmote  silt  loam 
Fairbanks  stlt  loam 
Fairbanks  silt  loam 
Mine  tailings 


Depth  to  bedrock 


Permafrost 

or 

table 

alluvial  gravel 

(m) 

(m) 

1 to  1.2 

More  than  1.8 

0.6  to  1.2 

1.2  to  1.8+ 

None  or 

0.3  to  1.8 

>4.6 

0 .3  to  1.8 

0.3  to  1 

More  than  1.8 

0.3  to  1 

1 to  1.8+ 

0.3  to  l 

More  than  1.8 

< 0.3 

0.3  to  0.6 

Variable 

More  than  1,8 

None 

More  than  1 

None 

0.3  to  0.6 

None 

0.6  to  1 

None 

1 to  1.8 

> 4.6 
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Figure  32.  Geological  cross  section  across  Little  Chena  Valley,  28  miles  along  the  Chena  Hot 
Springs  Road.  The  cross  section  is  based  on  drill  logs  provided  by  the  Alaska  District,  Corps  of 

Engineers. 

Correlating  the  soils  map  with  resistivity  demonstrates  that  the  decrease  in  resistivity  with 
frequency  is  associated  with  soils  having  organic  surface  layers,  or  partially  thawed  silt  covers; 
and  LF  and  BCB  values  in  excess  of  500  ohm-m  appear  to  be  associated  with  soils  without  organic 
surface  layers  or  without  completely  frozen  active  layers.  These  relationships  will  have  to  be 
examined  in  some  detail  in  future  field  studies.  The  mapped  Qsu  unit  in  the  northern  part  of  the 
study  area,  which  behaved  anomalously  at  VLF,  continues  to  behave  differently  at  LF  and  BCB 
compared  with  the  Qsu  unit  in  the  central  part  of  the  study  area. 

Clwu  Hot  Sprlap  Road 

This  small  study  area  is  in  the  flood  plain  where  the  variation  in  ground  conditions  often 
makes  interpretation  difficult.  The  sediment  types  found  vary  from  clean  gravels  to  silt  beds 
which  may  have  sections  that  are  continuously  frozen.  The  degree  of  interpretation  that  can  be 
made  from  airborne  resistivity  data  in  such  environments  depends  on  the  amount  of  ground  truth 
available,  and  at  this  site  unusually  good  subsurface  control  was  available.  The  Alaska  District, 
Corps  of  Engineers,  has  drilled  several  test  holes  across  the  Little  Chena  Valley  for  a preliminary 
investigation  for  a dam  site.  The  geological  cross  section  obtained  from  this  work  is  shown  in 
Figure  32,  and  the  location  of  the  cross  section  is  shown  on  the  resistivity  map  in  Figure  33. 

In  addition,  several  Wenner  resistivity  soundings  were  made  along  the  cross  sections. 


Figure  33.  The  contoured  resistivity  map  for  a sec- 
tion in  the  Little  Chen  a Valley,  28  miles  along  the 
Chena  Hot  Springs  Road.  The  dotted  line  indicates 
the  location  of  the  cross  section  in  Figure  32. 

The  contoured  resistivity  maps  of  the  three-frequency  airborne  survey  are  shown  in  Figure  33. 

The  L'?  and  BCB  survey  was  flown  across  the  valley  in  a N-S  direction  at  a 0.1-mile  line  spacing, 
while  the  VLF  survey  was  flown  at  a 0.8-mile  line  spacing  along  the  valley  in  a general  E-W  direc- 
tion. Figure  32  indicates  the  soil  and  geological  features  that  influence  resistivity  at  each  frequency. 

Figure  33  shows  that  at  the  southern  end  of  the  area  shallow,  conductive  bedrock  appears  to 
affect  the  VLF  data,  and  resistivities  less  than  250  ohm-m  are  encountered.  When  the  bedrock 
slopes  away,  the  VLF  resistivity  increases.  Farther  north  a silt  cover  blankets  the  gravel  and  its 
thickness  gradually  increases  to  the  north.  The  thawed  silt  blanket  in  the  central  part  of  the  area 
partially  screens  the  high  resistivity  in  the  gravel  at  VLF,  so  that  VLF  resistivities  are  limited  to 
from  500  to  1000  ohm-m.  The  areas  where  resistivities  at  VLF  are  in  excess  of  1000  ohm-m  prob- 
ably indicate  portions  of  the  section  where  the  silt  cover  is  entirely  or  partially  frozen. 

The  resistivities  at  LF  and  BCB  decrease  towards  the  north.  This  decrease  is  caused  by  a 
conductive  layer  at  the  surface;  the  conductor  is  probably  an  organic  surface  layer  and  a partially 
thawed  silt  cover. 

At  the  southern  end,  where  there  is  no  silt  cover,  the  resistivity  values  increase  with  frequency. 
The  skin  depth  at  I.F  and  BCB  is  insufficient  to  be  much  influenced  by  conductive  bedrock. 

Because  of  the  detailed  cross  section  available,  as  well  as  several  resistivity  measurements 
on  the  ground,  the  results  were  compared  with  model  studies.  Figure  34  shows  the  resistivity  pro- 
files along  the  cross  section  BB'  shown  in  Figure  33.  The  cross  section  here  was  divided  into 
3 sections,  and  in  each  section  at  least  one  Wenner  spread  was  made. 

A Wenner  spread  in  section  A is  shown  in  Figure  35.  A layer  of  gravel  d2  with  resistivities  in 
excess  of  1000  ohm-m  overlies  a conductive  stratum,  whicli  is  apparently  severely  decomposed  bed- 
rock, with  a resistivity p3  of  about  50  ohm-m.  The  very  high  resistivities  in  the  gravel  could  indi- 
cate frost  in  this  section  as  the  drill  logs  showed. 
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Chen  a Valley  and  the  geological  cross  section  derived  from  the  meas- 
ured resistivity  data  and  modeling  studies. 


Figure  35.  A Wenner  spread  in  the  Little  Chena  Valley,  in 
section  A of  Figure  34.  The  solid  line  represents  the  com- 
puted values  lor  the  ground  conditions  indicated. 


A Wenner  spread  in  section  B is  shown  in  Figure  20e.  A layer  of  seasonal  frost  in  the  silt 
cover  causes  relatively  high  resistivities  at  the  surface;  a subsequent  decrease  in  pB  is  caused  by 
a layer  of  unfrozen  silt  da;  the  silt  cover  in  turn  overlies  a layer  of  gravels.  Because  the  gravels 
are  of  lower  resistivity  here  than  in  section  A,  they  are  most  likely  umrozen. 

In  section  C,  part  of  the  silt  cover  may  be  frozen.  That  some  changes  occur  in  ground  condi- 
tions is  evident  from  a small  increase  in  both  the  VLF  and  LF  frequencies.  Because  of  an  ice 
layer  on  the  surface  in  this  section  in  April  1973,  the  Wenner  spreads  are  of  low  quality  because 
of  contact  resistance  problems. 
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Figure  36.  The  computed  values  of  apparent  resistivity  at  VLF  and  LF  in 
section  A,  Chen  a Hot  Springs  Road,  (or  the  E-phase  method  as  a function  of 
the  thickness  of  a highly  resistant  gravel  layer. 


The  situation  in  section  A is  modeled  in  Figure  36  at  VLF  and  LF;  pt  is  plotted  as  a function 
of  depth  to  bedrock.  In  accordance  with  the  measurement  of  the  Wenner  spreads,  the  gravel  dj  has 
been  given  a high  resistivity  of  104  ohm-m,  whereas  the  bedrock  is  given  a value  of  50  ohm-m. 
Figure  36  shows  that  a higher  value  for  pa  is  predicted  at  LF  than  at  VLF.  Figure  36  was  used  to 
predict  depth  to  bedrock  in  section  A,  and  the  predicted  depth  to  bedrock  is  plotted  in  Figure  34. 

The  drill  logs  show  that  a silt  blanket  begins  to  overlie  the  flood  plain  gravels  in  section  B. 
That  a new  geological  setting  occurs  is  also  evident  from  the  LF  resistivity  data  obtained  with 
airborne  equipment,  which  show  a reversal  in  trend.  In  Figure  22  the  effect  of  an  unfrozen  silt 
blanket  is  computed;  p&  is  given  as  a function  of  depth  to  bedrock  at  three  different  thicknesses 
of  the  silt  cover  at  VLF.  The  gravel  dg  has  been  assigned  a resistivity  of  1000  ohm-m  consistent 
with  the  data  of  the  Wenner  spread  in  Figure  20e.  The  depth  of  silt  refers  to  the  thickness  of  un- 
frozen silt.  In  Figure  23,  pg  is  plotted  at  LF  as  a function  of  the  thickness  of  the  silt  cover  dy 
The  apparent  resistivity  at  LF  becomes  insensitive  to  the  thickness  of  the  gravel  layer  and  the 
basement  resistivity.  From  Figure  23  the  thickness  of  the  silt  cover  can  be  estimated  and  from 
Figure  36  the  depth  to  bedrock.  Bedrock  at  a depth  greater  than  50  m cannot  be  distinguished. 

In  section  C part  of  the  silt  cover  is  frozen  and  part  is  unfrozen.  Both  VLF  and  LF  show  an 
anomaly  in  section  C.  It  is  difficult  to  model  this  complicated  situation  realistically.  When  the 
cross  section  shown  in  Figure  34,  obtained  from  the  resistivity  measurements  and  computer  model- 
ing, is  compared  with  the  cross  section  derived  from  drill  logs,  shown  in  Figure  32,  the  results  in 
sections  A and  B of  Figure  34  show  that  modeling  a simple  2-  or  3*layer  situation  yields  good 
agreement.  Such  a detailed  interpretation  can  only  be  made  if  good  ground  truth  data  are  available. 

Moose  Creek  Daa 

This  large  study  area  lies  almost  exclusively  on  the  alluvial  section  of  the  Chena  River.  The 
sediment  types  and  local  stratigraphy  are  inherently  variable,  including  interbedded  sands,  silts 
and  gravels.  Zones  in  the  finer  grained  material,  and  in  some  cases  the  sands  and  gravels,  may 
be  perennially  frozen. 


apparent  resistivity  contour  maps  at  VLF,  LF  and  BCB  in  the  Moose  Creek  Dam  area,  with  a preliminary  surficial  geological 
map  mainly  indicating  the  depth  of  silt  cover.  Clear  ••  250,  shaded  250-1000,  500,  solid  * 1000  ohm-m. 
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Figure  38.  The  computed  apparent  resistivity  at 
three  trequencies  for  two  layers  ot  silt  over  gravel 
as  a function  ot  the  thickness  of  silt  cover. 


The  area  was  flown  over  at  three  frequencies  and  represents  one  of  the  roost  difficult  areas  to 
interpret  because  of  the  mixture  of  sand,  silts  and  gravel  in  the  flood  plain.  The  contoured  resis* 
tivities  at  VLF.  LF  and  BCB  are  shown  in  Figure  37.  In  general,  resistivity  variation  in  this  area 
must  be  attributed  to  varying  thicknesses  of  silt  cover  over  poorly  graded  gravels  and  sands,  and 
the  presence  or  absence  of  perennially  frozen  sections  in  the  silt.  The  geological  map,  also  pre- 
sented in  Figure  37.  shows  the  preliminary  ground  data  provided  by  the  Corps  of  Engineers,  Alaska 
District.  The  solid  line  roughly  delineates  sections  of  varying  silt  cover  overlying  gravel: 

1.  In  area  A the  silt  cover  is  more  than  7 m thick  and  apparently  permafrost  is  frequently 

present  in  the  silt,  either  in  the  entire  depth  or  in  part  of  it. 

2.  In  area  B the  silt  cover  is  1.5  m to  5 m thick  and  only  sporadic  permafrost  is  present  in  the  silt. 

3.  In  area  C the  silt  cover  is  in  general  less  than  1.5  m thick  and  the  old  stream  channels  are 

occasionally  filled  with  silt. 

4.  In  area  D the  silt  cover  is  0 to  3 ra  deep. 

In  Figure  38  the  computed  apparent  resistivity  for  two  layers  of  silt  over  gravel  is  computed 
at  the  three  frequencies  used  in  the  survey.  The  apparent  resistivity  at  all  three  frequencies  de- 
creases with  increasing  silt  cover.  Very  little  difference  is  expected  between  data  at  LF  (275  kHz) 
and  data  at  BCB  (660  kHz).  For  this  type  of  layering,  the  resistivity  decreases  with  frequency. 

Thus,  from  the  modeling  curve  there  are  predicted:  a silt  cover  less  than  2 m thick  in  areas  where 
LF  and  BCB  values  are  in  excess  of  500  ohm-m,  and  a silt  cover  more  than  14  m thick  in  situations 
where  the  VLF  resistivity  (18.6  kHz)  is  less  than  250  ohm-im 

The  modeled  curves  refer  to  the  amount  of  thawed  silt,  and  frozen  layers  in  the  silt  will  prob- 
bably  be  indistinguishable  from  gravels.  Also,  bedrock  is  assumed  to  be  deep  enough  so  as  not  to 
affect  VLF  resistivities.  When  bedrock  comes  to  the  surface,  it  will  change  the  interpretation. 

In  comparing  the  resistivity  contour  map  with  boundaries  of  the  silt  cover,  several  areas  of  agree- 
ment are  observed: 

1.  Blocks  of  high  resistivities  (>  500  ohm-m)  are  observed  to  coincide  with  area  C and  the 
occurrence  is  consistent  at  all  three  frequencies.  Area  C has  a thin  silt  cover,  and  a thin  silt 
cover  is  also  predicted  from  the  modeling  curve.  The  consistency  with  frequency  in  this  area  is 
encouraging. 
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2.  A belt  of  low  resistivity  st  all  three  frequencies  is  observed  in  ares  A bat  also  extends  to 
the  north  of  it  in  ares  B.  From  the  modeling  curves,  silt  thicknesses  in  excess  of  IS  m are  expected 
because  of  the  low  VLF  resistivities  encountered.  A large  part  of  the  area  of  low  resistivity  coin* 
cides  with  area  A,  but  is  more  extensive  to  the  north  than  indicated  by  the  map  provided  by  the 
Alaska  District,  Corps  of  Engineers. 

3.  A block  of  higher  resistivities  is  again  encountered  at  LF  and  VLF  in  areas  B and  D.  The 
silt  cover  here  is  mapped  to  be  less  than  5 m thick. 

In  summary,  the  resistivity  data  in  this  area  are  difficult  to  interpret  because  of  the  complex 
and  random  geology.  However,  from  what  ground  truth  is  available,  a reasonable  agreement  between 
resistivity  contours  and  soil  types  is  observed. 


CONCLUSONB 

The  results  of  the  survey  in  the  four  study  areas  discussed  indicate  that  the  VLF  contour  maps 
show  the  most  apparent  correlation  between  surficial  geology  and  resistivity.  The  VLF  resistivity 
data  permit  a distinction  between  thawed  and  frozen  zones  on  fine-grained  sediments.  Both  in  Gold- 
stream  and  at  site  2 the  500  ohm-m  contour  line  often  closely  follows  the  mapped  geology  boundary  or 
thawed  and  frozen  units.  When  the  surficial  geology  is  complicated,  as  is  often  the  case  in  the 
flood  plain,  a correct  interpretation  of  resistivity  data  re  'ires  much  ground  truth.  However,  when 
good  ground  truth  data  are  available,  as  on  the  site  on  the  Cbena  Hot  Springs  Road,  much  informa- 
tion can  be  extracted  from  resistivity  data,  and  resistivity  data  provide  an  excellent  source  of  con- 
tinuous information  between  drill  holes. 

Although  the  LF  resistivity  data  are  useful  for  study  purposes,  the  additional  information  ob- 
tained from  the  data  would  have  to  be  weighed  against  the  extra  cost.  In  addition,  the  three  fre- 
quencies were  not  ideally  separated  in  this  survey.  Probably  the  ideal  frequencies  would  have  been 
VLF  and  a station  on  the  SLF  (70-140  kHz)  band.  An  SLF  station,  however,  was  not  available  in 
Fairbanks,  Alaska. 

Finally,  it  can  be  concluded  that  the  study  has  provided  an  adequate  test  of  multifrequency 
airborne  resistivity  mapping  in  the  discontinuous  permafrost  zone,  and  that  several  important  geo- 
technical objectives  can  be  derived  from  airborne  resistivity  mapping. 
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